Force and kinetic barriers to unzipping of the DNA double helix 
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■ A theory of the unzipping of double-stranded (ds) DNA is presented, and is compared to recent 
micromanipulation experiments. It is shown that the interactions which stabilize the double helix 
and the elastic rigidity of single strands (ss) simply determine the sequence dependent « 12 pN 
force threshold for DNA strand separation. Using a semi-microscopic model of the binding between 
nucleotide strands, we show that the greater rigidity of the strands when formed into dsDNA, 

■ relative to that of isolated strands, gives rise to a potential barrier to unzipping. The effects of this 
O ' barrier are derived analytically. The force to keep the extremities of the molecule at a fixed distance, 

the kinetic rates for strand unpairing at fixed applied force, and the rupture force as a function of 
loading rate are calculated. The dependence of the kinetics and of the rupture force on molecule 
length is also analyzed. 

■ Introduction: In cells, proteins apply forces to unzip and stretch DNA. These forces can be studied in single- 
^ molecule experiments (Fig. |^) 0-0], and are of biophysical as well as biological interest. Our focus here is primarily 

, on unzipping experiments where forces are applied across the double helix to adjacent 5' and 3' strands (Fig. |l]A) 
$3 |p|,^|j^-[lO|] . In experiments, the control parameters may be the force / itself, the distance between the last base pairs 
^ 2r, or the rate of force increase or 'loading rate' (Fig. |^). We discuss the results expected in all these situations. 
1^ We first use a thermodynamical equilibrium approach to show that the sequence-dependent force associated with 

unzipping of large DNAs, /„ ~ 12 piconewton (pN), can be simply deduced from the known free energy of DNA 
denaturation and the elasticity of single-stranded DNA. The unzipping experiments of Essevaz-Roulet et al (Fig. ^ 
101 and Rief et al (Fig. 0D) 0| are accurately described at this macroscopic level. 

Other experimentally observable aspects of unzipping can only be investigated using a more detailed description 
^ ' of base-pairing interactions. We therefore present a semi-microscopic model which accounts for hydrogen bonds 
^ ' and stacking interactions |jll|-|l^. We show that a free-energy potential barrier originates from the greater range of 
\^ • conformational fluctuation of DNA strands when isolated, relative to when they are bound together to form dsDNA 
\^ (Fig. 0). Our model can be investigated in detail and allows precise calculation of the effects of this barrier for the 
, initiation of unzipping and the kinetics of strand dissociation. 

' We compute the force necessary to keep apart the two extremities of the DNA molecules at some distance 2 r, as 
, well as the shape of the opening fork (Fig. ||). Due to the potential barrier, this force is much larger at small r (and 
' can reach some hundreds of pN) than the asymptotic value fu at large r. Analysis of unzipping in thermal equilibrium 
at the high level of precision possible in AFM experiments would allow unambiguous verification of this predicted 
force barrier. 

The barrier makes strand dissociation an activated process with dynamics that can be analyzed using nucleation 
^ . theory ||l^. Unzipping starts with a transition 'bubble' a few (< 4) bases long (Fig. 0). We calculate the free energy 

■ of this bubble, and determine how the dissociation rate depends on applied force and molecule length (Fig. ||). Results 
C , are compared to the experiments of Bonnet et al. and of Porschke ||l^ . 

O ■ Extending Evans' theory for the breaking of single bonds |l^ to the case of a one-dimensional polymer |l^, we 
^ ' then calculate the most probable rupture force when the DNA molecule is subjected to a force which increases at 
a constant 'loading rate' (Fig. 0). The dependence of the rupture force upon loading rate and molecular length 
could be quantitatively tested by AFM unzipping experiments; these results also shed light on AFM DNA-stretching 
experiments of Struntz et al. |^ and of Rief et al. Q . 
' Thermodynamic Description of Unzipping: An unstressed double helix is stabilized against spontaneously 
dissociating into its two strands by the interaction free energy per base pair, which from a thermodynamic perspective 
we may take to be some average amount g^. Although dependent on sequence, we may consider = —LAksT, the 
value determined from single molecule experiments on an AT rich sequence in A phage [0|, as a reference for the free 
energy difference between dsDNA and separated ssDNAs. Our emphasis is on an understanding of the free-energy 
balance in unzipping rather than to study inhomogenous sequence effects 

In the presence of applied torque F and unzipping force / (Fig. 0B) the free energy difference per base pair between 
unzipped and base-paired DNA strands is 
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Ag{T, /) = 5ssDNA - .gdsDNA = -ffo + Oq^ + 2 



(1) 



When Ag < 0, opening is thermodynamically favorable. The last two terms represent mechanical work done per 
base pair unzipping the double helix. In order, they are the work done by the torque {9o = 27r/10.5 is simply the 
change in strand winding angle during conversion of dsDNA to separated strands), and the stretching free energy of 
the unzipped single strands. 

The function gs{f) in (|l]) is the ssDNA stretching free energy per base at fixed force. The leading factor of 2 simply 
reflects the fact that two bases of ssDNA are created for each base-pair of dsDNA which is unzipped. The ssDNA 
elastic behavior is complicated by nucleotide- interaction effects [ p2| , but experimental force-extension curves for A 
phage ssDNA in 150 mM Na"*" are well described by a freely-jointed chain-like (FJCL) elastic response for forces > 1 
pN, with Kuhn length d = 15A |2Q,01. The corresponding free energy for forces up to w 20 pN is 



g^'^'^'if) = -ksT ^ log 



smhjd f/[kBT]) 
df/[kBT] 



(2) 



where the contour length per base pair is Iss = 27 /im/48.5 kb ~ 5.6 A. 

At zero apphed force, gs{0) = 0, thus ssDNA is stable when F < r„ = go/0o = —2AkBT, in good agreement with 
an experimental estimate of the twisting torque needed to denature an AT rich sequence in a A phage DNA |^ (the 
sign indicates a left-handed dsDNA- unwinding torque). In this case, the work done by the torque during opening is 
simply —go = lAksT per base pair. 

In the opposite case where torque is zero {i.e. for dsDNA with no constraint of its strand linking number), the 
critical unzipping force (at which Ag{f) = 0) is f^'"^^ = 11 pN. These results are in good agreement with the mid- 
range of unzipping forces encountered with experiments on inhomogeneous-sequence DNAs by Essevaz-Roulet et al. 
(Fig. |i]A: 12 pN threshold to start, then 10 to 15 pN during unzipping A-DNA 0). Data of Rief et al. (Fig. 0D; 20±3 
pN for poly(dG-dC), 9 ± 3 pN for poly(dA-dT) [Q) gives via (Q) denaturation free energies of 1.1 and 3.5 fcsT per 
AT and GC base pair respectively, in good agreement with thermodynamical data |2^. Finally, the projected length 
of one ssDNA nucleotide along the force direction at the unzipping transition is given as d^"^*^^ — 4A. Bockelman et 
al. used a similar theory to analyze unzipping force dependence on sequence jl]. 

The inset of Fig. ^ shows the curve in the torque-force plane on which Ag — 0, which is the 'phase boundary' 
separating dsDNA and unzipped ssDNAs. This boundary is predicted to have the shape /„ oc (F — TuY^"^ for small 

fu- 

For forces up to 15 pN, (^ is approximated to O.lSfcsT accuracy by the simple quadratic form 

5f (/) = . (3) 

where the ssDNA effective elastic constant is C = 0.12fcBT/A^. Using this form allows analytical solution for the 
unzipping force, — (C |5o|/2)^/^ — 12 pN 10 1^; at this force the projection of ssDNA monomer length along 



the force direction is = (2|(7o|/C)^^^ — 5A. This quadratic approximation is quantitatively nearly equivalent to 
the nonlinear model (^ at forces up to w 15 pN (e.g. note the accord between the torque-force 'phase boundaries' in 
Fig. 3, inset); this will be key to the continuum theory below. 

Semi-Microscopic Model of Strand Binding: Features at the nucleotide scale relevant to the onset of unzipping 
are ignored in models like (1). We therefore move to a model which uses the distances 2r(n) between corresponding 
?ith base pairs of the two strands, as degrees of freedom. The energy of the DNA strands is: 

H= rdn|im(r(n)) ("^ V + C/(r(n)) | (4) 



This model is similar to models previously used to describe thermal denaturation [ [12||l3| . The precise form of (Q) 
follows from our previously developed model for denaturation by torque Jl^ , p^ by integration over angular degrees 
of freedom, followed by continuum limit for the base index n. 

The first term in (^ describes interactions between neighboring bases along each strand, and must depend strongly 
on the inter-base half-distance r{n), since conformational fluctuations of the strands are highly quenched inside the 
double helix, relative to those along ssDNA. The strand rigidity is m (r) = E 6-''^''-^°^ + C, where i?o = 10 A is the 
double helix radius, and where 1/6 = 0.6 A is the separation at which the strand rigidity changes from its double-helix 
value E + C, to the much smaller ssDNA value C. We use the value C = 0.12 ksT/A^ from the previous section. The 
rigidity of the strands inside the double helix has been determined from Raman measurements of internal vibrations 
of dsDNA, toheE^ bSksT/A'^ H-H- 
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The second term in is a potential, acting between the two strands, made up of the hydrogen-bonding energy 
between corresponding bases, plus a torque energy: U{r) = t/jy(r) — TRi/r where i?i = 6 A We use the Morse 
potential form pl for the hydrogen bonding interaction UH{r) = L'[(e-'^(''^^) - 1)^ - 1], with D = bMksT 

and a = 6.3 A-^Fig. | inset). 

In thermal equilibrium, strand unpairing is described by the partition function Z = J I)r{Ti) g-H/ksT ^ 
be computed by use of a continuum transfer matrix technique along the n coordinate, leading to a Schrodinger-like 
equation: 



2m{r) dr"^ 



+ 



'4}{r) ^ g ilj{r) . (5) 



The free-energy potential V{r) = U {r) + [ksT / 2) \n[m{r) / m{oo)\ includes an entropic contribution due to the decrease 
in the rigidity m(r) with strand unbinding. This entropic potential (Fig. H) arises when going from the path integral to 
the Schrodinger equation [2l| ] with non constant mass | , and here it generates a large force barrier to the initiation 
of unzipping, which strongly affects the kinetics of strand separation. 

The lowest eigenvalue go is the equilibrium free energy, and the corresponding ipQ (r) describes thermal fluctuations 
of distance between the two strands. For 5 > the eigenvalue spectrum becomes continuous, corresponding to the 
appearance of completely separated ssDNAs at 5 = 0. When DNA is the thermodynamical favorable state 50 < 
represents the free energy per base pair of a long dsDNA, relative to separated ssDNAs. The fluctuations of r are 
confined to the Morse well, so 50 is well approximated if we take V{r) = y(i?o) + Unir) — Uh{Rq) and m{r) = m{Ro) 
for which (|^) is exactly soluble Inside the well, ^po{r) is the Morse ground state |2^]; outside the well it can 

be computed using the WKB approximation Using the parameters listed above, and at zero torque, the lowest 
eigenvalue is go = —lAksT, in accord with the corresponding number assumed in the previous section. Application 
of an unwinding torque gradually increases ga until at r„ = —2AkBT it becomes zero and the DNA unwinds, exactly 
as occurs in the thermodynamic model of the previous section (Fig. 3, inset). 

Force required to hold ssDNA ends at a given distance: This situation can be analyzed simply in terms of 
ijjoi^)! without further computation. This is because, for a semi- infinite dsDNA, the function '!/'o(^) corresponding to 
go is the probability distribution for the two ssDNA ends to fluctuate a distance 2r apart. Therefore the free energy 
associated with a fluctuation which separates the two ssDNA ends by a distance 2r (or equivalently the total work 
done separating the two ends to a distance 2r) is 

W{r) = -ksTlogMr) (6) 

up to an additive constant which is unimportant for our analysis. Note that, neglecting surface interactions, this is 
half the free energy associated to separate a distance 2r in the middle of a long dsDNA, since an interior 'bubble' is 
made of two 'forks'. 

The unzipping force that must be supplied to hold the two ssDNA ends a distance 2r (Fig. |^B) apart is thus just 
the derivative of (6), /(r) = dW {r) / d{2r) . This is in Fig. |[ which displays a large force barrier of w 270 pN as 
the strands are forced apart. The barrier peak occurs for a half-separation r — Ro = 0.5 A, and then decays to 
the long-molecule unzipping force /„ ~ 12 pN by r — i?o « 4 A. At large distances, the work done per base pair 
by unzipping is 2fudu — SksT, twice as much as the denaturation free energy —go] this is because the force must 
unzip the DNA, and extend the highly flexible ssDNAs. The peak force is large compared to the fluctuations in force 
associated with sequence [|l|. The force barrier will not be observable in large-scale unzipping experiments [Q, but 
should be observable in AFM studies. A stiff cantilever with roughly 0.1 A thermal noise should be used to measure 
the force barrier as a function of essentially flxed opening distance. 

In this fixed-distance experiment, one might also measure the shape of the opening 'fork', by determining the 
relation between opening distance and base position, n{r) (Fig. 0B). The most probable configuration n(r) satisfies 
the equation of motion associated with (^) , which expresses force balance along the chain, 

, , d^r 1 , , , / dr \ ^ , , , 

-(^);^ + 2"^WU) ^^^^^ 

We integrate (0) to obtain r(n), the shape of the opening 'fork'. Starting from the opening point where r — Ro and 
dr/dn = 0, n(r) = 2[E (r - Ro)/V'{Ra)]^/^. Far from the opening point n{r) = r/d^ ~ r/(5A) jjll. 

Unzipping kinetics at fixed force: Many experiments on short (10-100 bp) dsDNAs ( 17| and caption 

of Fig. |l|) probe the kinetics of strand separation. The equilibrium results discussed above are a starting point for 
a kinetic theory of unzipping based on nucleation theory 1 16 . The Schrodinger equation (|^) with the fluctuation 



corrected potential V{r) describes dsDNA and ssDNA as locally stable molecular states. The general problem faced in 
unzipping kinetics is the transition from an initially metastable state (dsDNA or ssDNA, depending on the force, see 
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below) to a final, stable (lower-free-energy) state. Strand dissociation requires the whole polymer chain to cross the free 
energy potential barrier of V{r) (Fig. 2), which makes the transition rate strongly length and force dependent. Using 
the effective potential V{r) (including the entropic barrier) corresponds to averaging over microscopic fluctuations of 
individual bases, restricting us to consider unzipping rates slow compared to those of these microscopic fluctuations 
(the experimentally relevant regime). 

The transition rate (equivalently the inverse lifetime) has the form familiar from transition-state theory fl^ : 

J. = e-«*/^-^ (8) 

It requires the activation free energy G* of a transition state, relative to the initial metastable state. The transition 
state is the saddle-point conflguration of the free energy, with one unstable direction leading monotonically down to 
the initial and final states, and is the dominant transition pathway [p^ . 

Our transition states are just the partially unzipped configurations r*{n) determined from (|^), for force-dependent 
boundary conditions consistent with the initial (rj) and final (rf) states. The activation free energy is 

G* ^ H[r*{Ti)] + {kBT/2) [ dnln[m[r* {n)]/m{oo)] ~ 2/r*(0) - 



dr^2m{r){V{r)~g,n) - 2f{rf-n) , (9) 

i.e. the free energy associated with (|^) minus the free energy of the metastable state (G,„ = Ng^n) from which the 
transitions occur. G* is in practice the free energy of the few-base-pair 'bubble' portion of the transition state which 
separates the unzipped and double-stranded regions. 

Finally, the rate prefactor vq is the linear growth rate of unstable perturbations around the saddle-point con- 
flguration We assume viscous dynamics with a friction coefiicient per base C = 67r?7i?o (water viscosity 
ry = 1 X 10^^ kg/(m • sec)). From a detailed calculation we find vq « -Da^/(4C) = 1 x 10^^ sec~^, essentially 
the ratio of the negative curvature of the Morse potential near the top of the well, to the friction coefficient. We 
now describe how to compute the dissociation rate and the nucleation bubble shape (i.e. boundary condition ri,rf) 
depending on the imzipping force. 

Kinetics of unzipping where ssDNA is stable (/ > /„): If a steady force is applied which is slightly bigger 
than the equilibrium unzipping threshold /„, then the initial dsDNA is metastable relative to separated strands: 
Ag{f) = —go — 2p/C < (0), and formula (||) is directly applicable (Fig. ^) to calculation of the dissociation rate 
The free energy of the initial metastable state is gm = go- The transition state in this case is a short ssDNA 
'nucleus' of n* bases at the open end (n = 0) of the dsDNA, with boundary conditions at the ends of fixed force / 
{dr/dn = 2//TO[r(n)]), and free energy equal to the metastable dsDNA value {V{r{n)) — 2p /m\r{n)] = go). This 
determines = r(n*) and rf = r{0). 

The size of the nucleation 'bubble' depends weakly on the force and is n*{f) ~ 4 bases, independent of the overall 
DNA molecule length. Therefore, the dissociation time i_ = l/v- (Fig. ^) is length- independent for forces above 
/,( = 12 pN. The decrease in dissociation time with forces > /„ is due to reduction of G* by the applied force. 
Beyond ~ 230 pN, the barrier in V{r) is completely overcome, and unzipping is immediate (t- = 1/vo)- Our 
computation addresses only the initial unzipping barrier-crossing event and does not include the time necessary to 
push the fork down the dsDNA, which would introduce a weak molecular-length dependence for forces > /„ in Fig. ||. 
Fork motion near fu may also be retarded by additional barriers introduced by inhomogeneous sequence. Kinetics of 
unzipping where ssDNA is metastable (/ < /„): Below the equilibrium unzipping force threshold, an infinitely 
long molecule is stable as a dsDNA. However, relatively short (< 20 bp) oligomers have a finite strand-dissociation 
time which can be accelerated by applied force. Here the final ssDNA state is metastable (although with a possibly 
long lifetime) relative to the initial dsDNA state. Therefore our calculation scheme applies not to unzipping (off-rate 
v-), but instead to annealing of ssDNAs (on-rate i/^, Fig. ||b). We compute the off-rate from the on-rate, using the 
equilibrium condition v-jvj^ = ^sU)I^bT ^ where N is the number of base-pairs in the molecule, and where Ag(/) 
is the force-dependent free energy difference per base pair between the paired and unpaired states . 

The calculation of the on-rate requires computation of a saddle point configuration with a nucleation bubble of 
dsDNA at the end of ssDNAs to which the force is applied (Fig. ^3). The relevant boundary conditions are analogous to 
those above, apart from the requirement that the ssDNAs be in close proximity but not base-paired (i.e. slightly outside 
the potential barrier of Fig. g) and the free energy of the metastable initial state equals g^. = 2 5s(/) = —"^P/C 
The dsDNA nucleation 'bubble' remains near 4 bp as the force / decreases from to zero, leading to an essentially 
TV-independent on-rate v+ (again the time necessary for the actual 'zipping up' is not included). Then, the N- 
dependence from the energy difference of paired and unpaired DNA results in a strong molecular- length dependence 
of V- as shown in Fig. || for forces less than fu ■ 
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At zero force, this becomes a calculation of dissociation time for free dsDNAs in solution with result t_ = 
-l^Q0.6Af-6.3 ggp_ rpj^g jj^gg^ pjg_ II sijo^g tjiig as a function of dsDNA length; a 10 bp dsDNA has a lifetime of 

roughly 1 sec; the exponential length-dependence results in a 30 bp DNA being stable for w 10^^ secw 30000 yr. Our 
estimate of t_ for N — 5, t- 0.3 msec, is in agreement with the results of Bonnet et al. (Fig. |i|e). The prediction 
for the nucleation bubble size, n* = 4, is close to the value n* — 3 measured by Porschke for a poly(A)poly(U) acid 

Rupture of dsDNA during gradual loading: The previous two sections discuss fixed-force experiments; an 
alternate experiment is to steadily increase force (/ = Xt where A is the 'loading rate' in pN/sec) and then to measure 
the force at which rupture occurs. Generally, rapid loading rates result in rupture at large forces. We are able to 
predict the most probable rupture force versus loading rate and molecular length, and our results display a rich range 
of possibilities (Fig. ^). 

Using the calculations of off-rates presented above, the probability distribution for rupture with force is 

Pif) - i._(7V, /) exp l^-i J^^ d/V_(A^, n] (10) 



For a number of molecular lengths. Fig. ^ shows the location of the peak of this distribution, the most probable 
rupture force that would be measured experimentally. For sufhciently slow loading rate, rupture occurs at zero force 
simply by thermal dissociation; for molecules > 20 bp thermal dissociation is practically unobservable. 

At some loading rate Ai, the peak in P(/) jumps to finite force and then increases, with each length of molecule 
following a different curve. In this regime, the forces at rupture are typically below the equilibrium unzipping threshold 
/„, since there is time for many thermal attempts at barrier-crossing to the metastable ssDNA state during loading. 
The length-dependence follows from the calculation of formation of metastable ssDNA discussed above, and has been 
qualitatively observed for 5'-5' pulling experiments (Fig. |l|C) ||^. 

At a higher loading rate A2 the peak of the rupture force distribution hits /„ which remains the most likely rupture 
force up to a loading rate A3. For A > A2, dissociation is occuring at forces large enough that the ssDNA final state 
is stable, resulting in no A^-dependence. Finally, beyond the very large loading rate A3 « 10^'^ pN/sec, the rupture 
force gradually increases simply because the molecule is unable to respond to the force before it becomes very large. 
A3 separates equilibrium and nonequilibrium time scales for very long sequences (N — > 00): the rupture force is 
independent on the loading rate and equal to /„ when A < A3, and increases above. 

Our framework treats nonequilibrium rupture of a one dimensional object, a development of previous theory p^ ] 
necessary for interpretation of unzipping experiments. In light of our detailed predictions, 3'-5' AFA/[ unzipping 
experiments (Fig. |l|A) should yield interesting results, and would be of help for an accurate determination of the 
free-energy potential V{r). 
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FIG. 1. Sketches of some experiments referred to in tlie text. All experiments are at room temperature and in physiological 
liquid buffers (PBS or Tris). Arrows symbolize the applied forces. A: Unzipping experiment of Essevaz-Roulet et al. [1]: 
the 3'-5' extremities of a A-phage DNA (49 kbp) are attached to a glass microscope slide (with translational velovity v = 40 
nm/sec) and a polystyrene bead connected to a glass microneedle (with stiffness k = 1.7 pN //im). The loading rate equals 
\ = kv = 0.06 pN/scc. When the force approaches 12 pN, the DNA starts to open. As unzipping proceeds, the distance 
between the two single strands extremities is controlled and the force varies between 10 and 15 pN depending on the sequence. 
B: parameters used in the theoretical description: force /, torque T and distance 2 r between the two single strands extremities. 
C: Stretching experiment of Strunz et al.[5]: a short ssDNA (10, 20 or 30 bp with about 60% GC content) is attached by 
one 5'-end to a surface, the complementary ssDNA is attached by the other 5'-cnd to an AFM tip. On approaching of the 
surface to the tip, a duplex may form that is loaded on retract until unbinding occurs. The distribution of the rupture forces 
is obtained for loading rates ranging from 16 to 4000 pN/sec. D: Stretching and unzipping experiment of Rief et al.[4]: DNA 
of poly(dA-dT) (5100 bp) or poly (dG-dC) (1260 bp) are attached between a gold surface and an AFM tip, and stretched. 
Through a melting transition, single DNA strands arc prepared; these strands upon relaxation reanneal into hairpins as a 
result of their self-complementary sequences. The forces of unzipping of these hairpins are 20 ± 3 pN for poly (dG-dC) and 
9 ± 3 pN for poly(dA-dT). E: Dissociation experiment of Bonnet et al. [7]: The rate of unzipping, V-, and closing, of a 
5 bp DNA hairpin (CCCAA-TTGGG) is investigated by fluorescence energy transfer and correlation spectroscopy techniques. 
The hairpin is closed by a loop of 12 to 21 Thymine (T) or Adenine (A). The characteristic time of opening t- = l/v- is found 
to be largely independent of the loop length, and equal to t- ~ 0.5 msec. 
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FIG. 2. Base pair potentials in unit of fesT, as a function of the base radius r (in A), without (inset) and with (main picture) 
entropic contributions. Inset: Morse potential J7(r) accounting for the hydrogen bond interaction. Main picture: total potential 
V{r) for zero torque. Once entropic contributions are considered, small r values are less favorable and a barrier appears. The 
free-energy gdsDNA = go of the dsDNA is lower than the single strand free-energy Qssdna = 0. Note the difference of scales on 
the horizontal axis between the two figures. 
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Cocco. Fig 3. Desired size: 1 column, height 8.5 cm. 



FIG. 3. Force /(r) (in pN) to be exerted on the DNA to keep extremities at a distance 2r apart (in A). The peak force 
/ ~ 270 pN, reached at r ~ 10.5 A, is much larger than the asymptotic value ~ 12 pN, equal to the equilibrium force fu (at 
zero torque, and in the Gaussian approximation) for unzipping a large portion of the molecule. Inset: phase diagram, in the 
plane of torque F (fcsT) and of force / (pN). The lines shows the critical unzipping force /„ as a function of F with formula 
(2) (full line) and formula (3) (dashed line) for the stretching free energy of the single strand. Below the line, dsDNA is the 
stable thermodynamical configuration while for forces larger than /u(F), denaturation takes place. /„ vanishes at the critical 
torque Fu ~ —2.4 ksT, and is equal to 11 pN (full line) and 12 pN (dashed line) at zero torque. 
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Cocco. Fig 4. Desired size: 2 columns, height 6.5 cm. 



FIG. 4. Transitions states involved in the theoretical calculation of the kinetic rates. A - Unzipping: opening of dsDNA 
is favorable at forces f > fu, and the unzipping rate v- is calculated directly. The nucleation bubble is of (~ 4) base pairs, 
weakly depending on the force. B - Annealing: when f < fu, dsDNA is thermodynamically stable; the dissociation rate v- is 

obtained indirectly through the calculation of the annealing rate i/+ of the motastablo ssDNA, i/_ = i/^e~'^'^^^^^^'^^'^ , where 
Ag{f) > is the excess of free-energy per bp of ssDNA with respect to dsDNA. The nucleation bubble is of (~ 2) base pairs. 
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FIG. 5. Time of dissociation f_ (in sec) as a function of the force / (in pN). Three regimes can be distinguished. For 

/ < /„ = 12 pN, the dissociation times depend on the length A'^ of the sequence (A'^ = 10, 20, 30 bp from bottom to top). For 
fu<f<fb = 230 pN, the dissociation time is length independent and decreases, as the energetic barrier to overcome, lowers. 
For f > fb, no barrier is left and dissociation is immediate. The slope of the logarithm of t- near /„ is d\ogt- /df = —8 
^ {f > fu), —2duN + 31 A with du = SA (/ < /„). Inset: Time of thermal dissociation t- (for zero force) as a function of the 
number of base pairs N . 
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FIG. 6. Rupture force (pN) as a function of the loading rate A (pN/scc) for five different molecule lengths A'^=10, 20, 30, 
50 and 100. Arrows indicate the different critical loading rates for A'^ — 10. Below Ai (=10 " * for N — 10), rupture occurs at 
essentially zero force through thermal dissociation. For loading rates ranging from Ai up to A2 (=10'*'® for A'' = 10), the rupture 
force is finite, and thermal tunneling is responsible for the strong dependence on N, until the force reaches the equilibrium value 
fu = 12 pN. For larger loading rates, the rupture force is length independent. It increases again as A > A3 = 10^'^ pN/sec, 
since the molecule is unable to respond to the force before it becomes very large. 



12 



